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I. INTRODUCTION 



In order to extract fundamental quantities from systems containing heavy quarks, such as 
the weak mixing angles Vf, c and V^, the strong interaction contributions to the observable of 
interest must be determined. This is because the standard model of electroweak interactions 
is constructed in terms of quarks, gluons and leptons, while we are only able to study the 
quarks via the properties of the hadrons within which they are confined. Significant progress 
was made in eliminating the strong interaction uncertainties in some matrix elements in 
hadrons containing one or more heavy quarks in the late eighties with the discovery of 
the heavy quark symmetries |TJ], and the formulation of heavy quark effective theory 0. 
Ultimately, such strong interaction effects will be directly computable from numerical lattice 
QCD simulations, however, at this point in time, and in the foreseeable future, lattice 
simulations will not be performed with the physical values of the light-quark masses, m q . 
In order to make a connection between the lattice simulations and nature, an extrapolation 
from the lattice quark masses to those of nature is required. One can most efficiently 
perform the extrapolation for many observables using simulations of partially quenched 
QCD (PQQCD) matched to partially quenched chiral perturbation theory (PQ^PT) 0-0. 
Some observables have been analyzed with PQ^PT in the light pseudo-scalar meson sector, 
such as decay constants and masses |J. Further, the decay constants and B-parameters for 
mesons containing a single heavy quark have been computed . In this work we present the 
one-loop analysis in PQxPT of the Isgur-Wise function for B — > decays, of the matrix 
elements of isovector twist-2 operators that are directly related to moments of the parton 
distribution functions, and of the radiative decays D* — > Z>y. 



II. PQQCD 

We will consider a theory constructed from three valence- quarks, u, d, s, three-sea quarks 
j, I, r and three bosonic-quarks u, d, s. The masses of the bosonic-quarks are equal to those 
of the valence- quarks, m, m, m s respectively, where we have chosen to work in the isospin 
limit. The sea-quarks are taken to be degenerate with mass mj. The valence-, sea-, and 
bosonic-quarks are combined into column vectors 

Q L = (u 7 d,sJ,l 7 r,u,d,f) L , Q R = (u,d,8,j,l,r,u,d,ij T R , (1) 

which transform in the fundamental representation of SU(6\3)L t R respectively. The ground 
floor of Ql transforms as a (6, 1) of SU(6) q L <8> SU{3)qL while the first floor transforms as 
(1,3). The right handed field Qr transforms analogously. The PQQCD Lagrange density 
is invariant under SU(6\3)l ® SU(6\3)r chiral transformations, and in direct analogy with 
QCD a chiral Lagrangian can be constructed to describe the low-energy dynamics of the low 
lying hadrons. 

The low-energy QCD dynamics of hadrons containing a single heavy quark are described 
by a chiral lagrangian that has the heavy quark spin symmetry and flavor symmetry man- 
ifest [10,11]. The extension of the heavy quark chiral Lagrange density to describe heavy 



quark systems in PQQCD can be found in Ref. HJ. The dynamics of the pseudo-scalar 
mesons are described at leading order by a Lagrange density of the form, 
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where M and M are matrices containing bosonic mesons while x an d X^ are matrices con- 
taining fermionic mesons, with 
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where the upper 3x3 block of M is the usual octet of pseudo-scalar mesons while the 
remaining entries correspond to mesons formed with the sea-quarks. The singlet field is 
defined to be $o — str ( $ ) / "\/2, and its mass mo can be taken to be of order the scale of 
chiral symmetry breaking, mo — > A x ||. The super mass matrix, rriQ, is 



rriQ = diag (m, m, m s , rrij, rrij, rrij, m, m, m s 



(5) 



where we will work in the limit of exact isospin symmetry. The convention we use corresponds 
to / ~ 132 MeV. 

The B-mesons with quantum numbers of bQ form a nine- component vector, 



B = (B u , B d , B s , Bj, Bi, B r , Ba, B s 



(6) 



and heavy quark spin symmetry (for a comprehensive review see Ref. ||16|| ) is implemented 
by combining the annihilation operators for the B and B* mesons together into the field H v 



H^-il + ^i^Bl + i^B 



(7) 



where H v — > H V V^ under SU(3) flavor transformations and H v — > SbH v under heavy quark 
spin transformations. The four-vector is the four-velocity of the heavy meson. The low- 
momentum strong interaction dynamics of the heavy mesons are described by a Lagrange 
density of the form M , 
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where 
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and where the light-meson fields are 
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The Tr [ ] 's in eq. ([8]) correspond to traces over Dirac indices, and implicit in eq. (||) are 
str [ ]'s over the flavor indices. Couplings to $0 H, such as str [ ), have not been included 
as the $0 will be integrated out of the theory, and such contributions will be included via 
higher dimension local operators. The axial coupling constant g n has been constrained to be 
g n ~ 0.56 JT2| (or g w = 0.24) by the radiative decays D* — > and recent quenched lattice 

simulations yield g n ~ 0.42 p^JT5[] . The low-momentum strong interaction dynamics of 
mesons are described in an analogous way, with heavy quark flavor symmetry dictating 
that the axial coupling, g n , for the 5's and .D's is the same. Further, the dynamics of the 
anti- heavy mesons can be described in analogous way [ffl]. 



III. DECAY CONSTANTS 



The heavy meson decay constants, such as Jb u , have been analyzed at the one-loop level 
in PQxPT by Sharpe and Zhang Their analysis provides a very clear demonstration of 
the implementation of PQQCD and PQxPT, and it is worth reviewing their results. 

The decay constants of the heavy mesons are defined via the matrix element 



(0|g 7 M (l - lB)b\B q (v)) = -i f Bq m B ^ 



(11) 



In the heavy meson effective field theory, this matrix element receives contributions from 
both tree- level operators and from diagrams involving light-meson loops. Up to O (m q ), the 
local operators that contribute are 
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where n a is a column vector that picks out the light quark associated with the de- 
sired current, A is the renormalization scale, and M.+ is related to the mass matrix via 
M.± = \ (^itlqC, ± ^ffig^j . At leading order one finds that fs q = k/ y/iriB q - At O (m q ) in 
PQxPT the decay constants are M 
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1 A more detailed study @ gives g w ~ 0.76 ± 0.2 or g w ~ 0.27 ± 0.07. 
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where m qq , rnjj and rrij q are the masses of light mesons with quark composition qq, jj 
and jq respectively. The Cl R \A) are renormalized constants where both divergences and 
terms that are analytic in m q have been absorbed. The expression in eq. ( |T3| ) holds for all 
q = u,d, s,j, I, r, and agrees with the expression in Ref. J4| when there are three sea quarks, 
Nf = 3. Lattice simulations of f Bu d s and f Bj at various values of the sea-quark masses, mj, 
and various values of the valence masses will allow for the determination of the quantities 
k, C3 (A) and the combination C± (A) + C% (A). As these counterterms have the same 
value in QCD and PQQCD, their values can be inserted into the QCD expressions for the 
decay constants [[|.[Iy]] 
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IV. CHIRAL 1/Mg CORRECTIONS TO i? — > DW AT ZERO RECOIL 

The matrix elements for the semileptonic decays B — > Dez/ and B — > D*ez/ at zero-recoil 
are normalized to unity in the heavy quark limit Further, there are no corrections at 



zero-recoil of the form l/m c ,b by Luke's theorem ||18|| , and the leading corrections in the heavy 
quark expansion enter at l/m 2 . In addition to contributions from perturbative insertions 
of local operators in the heavy quark Lagrange density, there are also contributions that 
receive infrared enhancements from long-distance strong-interactions processes. The largest 
contribution of this type is from the one-loop diagrams shown in fig. p] in which the D* — D 
mass-splitting is explicitly retained in the propagator fl9[| . Recently, quenched lattice 
QCD simulations of the B — > form factors have been performed |20],^ll but the light- 
quark mass dependence of the results has yet to be ascertained.. 
The hadronic matrix elements for B — > D^eu are 

(D(t/)\&fb\B(v)) = ^m B m D [h + {w){v + v'Y + h-(w)(v - vJ L } 
(D*{v', e)\c^ l5 b\B{v)) = y/m B m D * [ -ih Al {w){w + l)e*" + ih Aa {w)(v ■ e>" 

+ih A3 (w)(v e*)v^ ] 

(D*(v', e)\crfb\B{v)) = ^m B m D * h v {w)e^elv' a v p , (15) 

where is the four-velocity of the initial-state S-meson, v'^ is the four-velocity of the final- 
state D^*> -meson, w = v ■ v' and h±, hy and h Al2i are the six independent form factors 
that contribute. In the heavy quark limit h_{w) = h A2 {w) = and h + {w) = h v {w) = 
h Al3 (w) = where £(w) is the universal function for the B^ — > decays known 

as the Isgur-Wise function. The matrix element is reproduced in the heavy meson chiral 
lagrangian by an operator of the form 
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FIG. 1. One-loop diagrams that contribute to the matrix elements for B —* D^*>. The solid 
lines denote mesons containing a heavy quark, while the dashed lines denote light mesons. The 
solid square corresponds to an insertion of the weak interaction vertex, while the small solid circles 
denote a strong interaction vertex oc g^. Diagram (a) is a vertex correction, while diagrams (b) 
and (c) correspond to wavefunction renormalization. 



! (l-75)^ 6) 



(16) 



where £,{w) is the Isgur-Wise function, and we have retained the flavor labels on the H fields. 

The lowest order heavy-quark-spin-symmetry breaking operator in the heavy quark chiral 
Lagrangian is dimension three due to the chromomagnetic moment interactions of the charm 
quark with the light degrees of freedom, ~ c a^G^ cjm c (we neglect the contribution from 
the b-quark chromomagnetic interactions, i.e. rrib 3> m c ). This interaction gives rise to a 
mass-splitting between the D* and D at order l/m c in the heavy meson chiral Lagrangian, 
which after performing the Dirac trace, is 



(17) 



As this is the lowest dimension operator that breaks the heavy quark symmetries, it provides 
the dominant long-distance contribution to the 1/m 2 corrections to the matrix elements in 

eq. ( |ToD |H| . Randall and Wise JHJ found a contribution of the form log (m 2 /A 2 ) + ... 

at w = 1 from the pion loop-diagrams shown in fig. |l|. In order to exactly compensate the 
renormalization scale dependence, A, the 1/m 2 operators that are independent of m q will 
combine together to give a m g -independent contribution X + (A) at w — 1 for the B — > D 
decay and Xa{A) at w = 1 for the B — > D* decay. The form factors at zero-recoil are 
modified at O (1/m 2 ) by the one-loop diagrams in fig. || to become 
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where we have used the shorthand F ab = F(m ab , A^ c \ A), F ab = F(m ab , — A^, A), R ab = 
R(m ab , A^ c \ A) and R ab = R(m ab , — A^, A). The expressions in eq. flTSp are true for 
g = it, d, s. The loop functions F and R are 
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It is interesting to note that in the A^ c ^ <C rn qq limit, there is a contribution to the h\ Bq ^ P ®(l) 
of the form (m 2 q — m 2 j)/m 2 qq . This is analogous to the enhanced chiral logarithms found by 
Sharpe || that contribute to the light meson masses and decay constants. Once lattice sim- 
ulations have determined X + (A) and Xa{A), they can be inserted into the QCD expressions 
for h + (l) and h Al (l) for B u ^ s to obtain the zero-recoil matrix elements in QCD, 
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where we have used Fy = F(my, A^ c \ A) and Fy = F(my, — A^ c \ A). Keeping only the 
contribution from pion-loops and counterterms we recover the results of Ref. 



V. CHIRAL CORRECTIONS TO THE B -> D ISGUR-WISE FUNCTION 



While the Isgur-Wise function £(w) in eq. (|TB[) is normalized to unity at zero recoil, £(1) = 
1, away from zero recoil its value depends upon the light quark masses. In QCD the leading 
non-analytic contributions to the Isgur-Wise function for B u ^ decays, £b u d (w), and the Isgur- 
Wise function for B s decays, £b s (u>), of the form m q \ogm q have been computed p2| , p3| . In 
addition, the contributions of this form arising in quenched QCD (QQCD) have also been 
computed pi}] . In the heavy quark limit, matrix elements of the operator C7 M (1 — 75)6 



between B and D^*> states are reproduced, up to 0(m g ), by 
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where we use notation similar to Ref. [^3] for the higher order counterterms, 771,2, and it 
is understood that flavor indices are super-traced over. In addition heavy quark symmetry 
requires that 771,2(1, A) = The one-loop diagrams shown in fig. |T] and the contributions 
from eq. ( pl|) give Isgur-Wise functions in QCD at O (m q ) 
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where the w-dependent function is [[22| , |23 
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7712 are renormalized constants with 77^(1, A) = 0. 

In PQxPT, the Isgur-Wise functions for mesons composed of a b-quark and an anti- 
valence quark, bu, bd, or bs, at one-loop order and in the limit of isospin symmetry, are 
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where q = u,d, s. The Isgur-Wise function for mesons comprised of a b-quark and an anti-sea 
quark, j, is 
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As required, the expression in eq. (|22| ) reduces to the expression in eq. (|25| ) when 
m^ = rrijj and m u ^, s = nij. By determining and ( w ) i n lattice simulations and 

fitting the quantities £ (w) and 771,2 (w,A), one can then use the QCD expressions for £ Bu d s (w) 
in eq. (^) to perform the extrapolation in m q . 



VI. RADIATIVE TRANSITIONS: D* -> L> 7 

The widths of the radiative decays D* — > are relatively well measured. These decays 
are quite interesting as they have allowed for a determination of the axial coupling constant 



g n , in combination with the leading order magnetic transition moment counterterm [12,13 
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The value for g n that has been extracted from these decays is significantly smaller |T2|JT3 
than the value expected in the naive quark model. 

The width for radiative decays is given in terms of a total radiative matrix element fi a , 



T(D:^D al ) = -\^ a \ 2 \k 1 \ 3 , (26) 

where the subscript "a" denotes the flavor of the meson. At leading order there are 
contributions that are independent of the light quark masses arising from the coupling of 
the photon to the light degrees of freedom, and to the heavy quark. At next to leading 
order there is a contribution of the form ^ffn^ resulting from the one-loop diagrams shown 
in fig. In the heavy quark limit, the transition moments in QCD are 

D* D 
« f 

\ / 




FIG. 2. One-loop diagrams that contribute to the radiative transitions D* — > D*y. The solid 
lines denote mesons containing a heavy quark, the dashed lines denote light mesons and the small 
solid-circle denotes a strong interaction oc g n . 
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where (3 is the m g -independent counterterm contribution from the light degrees of freedom. 
In PQQCD the calculations are analogous to those of QCD, and one finds that 

,.PQ _ 2 2 gl m ju 

^ " 3m c 3 P 2vr/ 2 

pq _ 2 1 glm ju 

^ ~ 3m c 3^ + 4vr/ 2 

^ s ~ 3m c 3 P + 4tt/2 

where the coefficient (3 is the same in QCD and PQQCD. A fit to lattice calculations of 
these matrix elements will provide a determination of (3, and consequently an estimate of 
the radiative matrix elements for the physical values of m q . 
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VII. MATRIX ELEMENTS OF ISOVECTOR TWIST-2 OPERATORS 



Recently, it has been realized that the chiral corrections to the matrix elements of isovec- 
tor twist-2 operators in the nucleon and octet-baryons can be computed in chiral pertur- 



bation theory |25|-p9|. Such corrections systematically incorporate the long-distance strong 
interaction contributions to the moments of the parton distributions. While it is not possi- 
ble to perform deep-inelastic-scattering (DIS) from hadrons containing heavy quarks, it is 
possible to consider the parton distributions of heavy hadrons from a theoretical standpoint, 
and further it is likely that lattice studies of these distributions will be significantly easier 
than analogous studies in nucleons. 

In QCD the twist-2 isovector operators have the form 
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where the {...} denotes symmetrization on all Lorentz indices, and A a are SU(3) Gell-Mann 
matrices. In PQQCD the operator has an analogous form, 
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where A is a super Gell-Mann matrix. As an example, the A matrix has entries 
diag(l, —1, 0, 1, —1, 0, 1, —1, 0). Forward matrix elements of Offl' 3 „ between heavy meson 
states are reproduced up to 0(m g ) by forward matrix elements of 
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where A 3 



\ (<^ 3 £ + ^ 3 ^)> an d it is understood that flavor indices are super-traced 
over. At O (m q ) in the chiral expansion the forward matrix elements receive contributions 
from the operators in eq. (|31|) and the one-loop diagrams shown in fig. [3|, 
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for n = odd and n > 1. In the loop calculations we have only retained the non-analytic 
contributions, and the are renormalized counterterms. For n = 1 the matrix element 

is absolutely normalized since O^' 3 is the isovector charge operator. Matrix elements in the 
other meson states can be found straightforwardly from those in eq. (|32f) . From PQQCD 
lattice simulations one envisages that the constants t^ N \ Sg (A) and the combination 
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FIG. 3. One-loop diagrams that contribute to the matrix elements of the isovector twist- 2 oper- 
ators in mesons containing a heavy quark. The solid lines denote mesons containing a heavy quark, 
while the dashed lines denote light mesons. The crossed- circle denotes an insertion of C^/l" — Mn 
while the small solid-circle denotes a strong interaction oc g w . 



expressions in QCD to determine the moments of combinations of parton distributions in 
heavy mesons. In QCD I find the matrix element between B u states to be 



(B |0(")' 3 \B ) 



t 



(n) 



:i+3<£) 

167T 2 / 2 



2 ,„ m l , ™2 ^rn 2 K 



2m$ log + m^log ^ 



+ ( s («)W ( A ) + 4«)H ( A ) + 2s f)W (A) 



. (R)(n) 

m + S3 



(A) m s 



(33) 



VIII. CONCLUSIONS 



In this work I have computed several observables in mesons containing heavy quarks at 
one-loop level in partially quenched heavy quark chiral perturbation theory and in ordinary 
heavy quark chiral perturbation theory. In order to determine these observables from nu- 
merical lattice QCD simulations, an extrapolation in the quark masses from those that can 
be simulated on the lattice down to their physical values is required. The recent progress in 
understanding partially quenched QCD has provided an efficient pathway to perform such 
extrapolations. The counterterms that contribute to partially quenched observables are 
the same as those that contribute to the analogous QCD observables, and hence partially 
quenched simulations along with the light-quark mass dependence about the chiral limit are 
sufficient to recover the QCD prediction to some order in the chiral expansion. 
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